Emerging probing perspective of two-dimensional materials physics: terahertz emission spectroscopy

Terahertz (THz) emission spectroscopy (TES) has emerged as a highly effective and versatile technique for investigating the photoelectric properties of diverse materials and nonlinear physical processes in the past few decades. Concurrently, research on two-dimensional (2D) materials has experienced substantial growth due to their atomically thin structures, exceptional mechanical and optoelectronic properties, and the potential for applications in flexible electronics, sensing, and nanoelectronics. Specifically, these materials offer advantages such as tunable bandgap, high carrier mobility, wideband optical absorption, and relatively short carrier lifetime. By applying TES to investigate the 2D materials, their interfaces and heterostructures, rich information about the interplay among photons, charges, phonons and spins can be unfolded, which provides fundamental understanding for future applications. Thus it is timely to review the nonlinear processes underlying THz emission in 2D materials including optical rectification, photon-drag, high-order harmonic generation and spin-to-charge conversion, showcasing the rich diversity of the TES employed to unravel the complex nature of these materials. Typical applications based on THz emissions, such as THz lasers, ultrafast imaging and biosensors, are also discussed. Step further, we analyzed the unique advantages of spintronic terahertz emitters and the future technological advancements in the development of new THz generation mechanisms leading to advanced THz sources characterized by wide bandwidth, high power and integration, suitable for industrial and commercial applications. The continuous advancement and integration of TES with the study of 2D materials and heterostructures promise to revolutionize research in different areas, including basic materials physics, novel optoelectronic devices, and chips for post-Moore’s era.


Introduction
Within electromagnetic spectrum, the THz region, historically spanning frequencies from 0.1 THz to 10 THz, has been referred to as the "terahertz gap" between microwaves and far-infrared radiation, posing challenges in generating and detecting radiation within this range (Fig. 1a) 1 .However, recent advancements in technologies such as optical rectification using femtosecond laser pulses, quantum cascade lasers (QCLs) 2 , tunable THz waveguide [3][4][5] and free electron lasers 6 have significantly expanded access to the THz part of the electromagnetic spectrum, overcoming previous limitations and creating new opportunities for research and applications in fields including THz spectroscopy, material science, communication, agriculture, environment, medical imaging, security screening and food [7][8][9][10][11][12][13][14] .Particularly in the field of material science, THz spectroscopy has been successfully employed to characterize different physical properties of materials, such as complex dielectric constant, refractive index, conductivities and so on.The unique frequency range of THz spectroscopy makes it well-suited for investigating intraband electrical transport and lowenergy excitations in emerging materials.Its contact-free nature as a free-space optical technique eliminates the need for electrical contacts and minimizes material contamination or modification.Additionally, the ultrashort duration of sub-picosecond THz pulses enables timeresolved studies of ultrafast processes using techniques like TES, optical pump-THz probe (OPTP) spectroscopy or THz pump-THz probe (TPTP) spectroscopy, among which TES has emerged as a powerful tool for investigating materials' electronic and optical properties.By harnessing the unique properties of emitted THz waves, researchers can delve into the intricate interplay between photons, charges, phonons, and spins (Fig. 1b) 15 , enabling the exploration of complex physical phenomena and intrinsic properties in materials.These advancements in THz technology not only bridge the historical 'terahertz gap' but also open new avenues for scientific research and practical applications, laying the groundwork for the next sections to discuss the specific uses and implications of these technologies in various fields.
Meanwhile, 2D materials have garnered significant attention for their extraordinary electronic, optical, and mechanical properties [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30] .Graphene, transition metal dichalcogenides (TMDs), black phosphorus, quantum wells, and 2D electron gas, along with additional examples such as hexagonal boron nitride (hBN) and layered perovskite, stand out as prominent examples of 2D materials [31][32][33][34][35][36][37][38][39] .Due to their atomically thin structure, 2D materials exhibit high surface-to-volume ratios, enabling enhanced surface interactions and sensitivity to external stimuli, which open up new possibilities for sensing, catalysis, and energy applications.Furthermore, the integration and stacking of 2D materials in van der Waals (vdW) heterostructures have significantly enhanced artificial lattice structures, leading to the formation of novel material systems that exhibit unique and unprecedented physical properties (Fig. 1c) 40 .By stacking different 2D materials and controlling their stacking angles, one can engineer and tailor the electronic and optical properties that allow the discovery of intriguing phenomena such as  15 .c Reprinted with permission from ref. 40 .d Reprinted with permission from ref. 52 .e Reprinted with permission from ref. 53 .f Reprinted with permission from ref. 54 .g Reprinted with permission from ref. 56 .Fig h, reprinted with permission from ref. 57 superconductivity, interfacial charge transfer, interlayer exciton formation and so on.From an industrial perspective, notable progress has been achieved in the fabrication [41][42][43][44] , device development [45][46][47] , and industrial applications of 2D materials [48][49][50] .These advancements position 2D materials to have a pivotal role in the post-Moore's law era of the chip industry 51 .Consequently, the development and integration of these materials into devices are essential for technological advancement and shaping the future of the semiconductor industry, surpassing the constraints of conventional silicon-based chips.To fully understand 2D material physics and use the full capabilities of those materials requires sophisticated characterization techniques, while the TES has been shown to be an excellent optical approach to studying 2D materials (Fig. 1d) 52 with non-intrusive probing of systems.In principle, TES measures the emitted THz time domain electric field that contains both amplitude and phase information, which can be used to fully detect the structure and electrical properties of materials.For example, pioneering research has been conducted using TES to study 2D materials such as MoS 2 monolayer (Fig. 1e) 53 , which reveals the unique play of surface depletion field-induced photocurrent surge in addition to optical rectification.In 2D heterostructures, interfacial effects dominate the nonlinear processes, such as active THz surface emission in Graphene/SiO 2 /Si heterostructures (Fig. 1f) 54 , THz wave modulation via metasurfaces in MoS 2 /Si heterostructures 55 , and interfacial current triggered THz emission in WS 2 /MoS 2 heterostructures (Fig. 1g) 56 .Besides, the complex physical phenomena and mechanisms such as Rashba splitting (Fig. 1h) 57 are also been deeply studied.Those different materials and heterostructures serve as fascinating platforms to hold different THz emission mechanisms, thereby enabling the exploration of complex physical phenomena and intrinsic material properties.Firstly, 2D materials provide high mobility and ultrafast carrier dynamics, enabling efficient generation and detection of THz waves.Secondly, the strong light-matter interaction in 2D materials allows for enhanced THz emission, making the spectroscopy more sensitive and accurate.Additionally, the flexibility and tunable electronic properties of 2D materials allow for the development of compact and versatile THz devices.This can lead to advancements in a wide range of applications, from security scanning and medical imaging to wireless communication and material characterization.The integration of 2D materials into THz spectroscopy thus opens up new avenues for exploring and exploiting the THz spectrum in both fundamental research and practical applications.
Hereby we comprehensively review TES and its application in studying 2D materials and heterostructures.Particularly we focus on the diverse nonlinear processes underlying THz emission in 2D materials and interfaces, such as optical rectification, photon-drag, high-order harmonic THz generation, and spin-to-charge conversion.By utilizing the capabilities of TES, researchers gain deeper insights into the electronic and optical properties of 2D materials, opening the way for progress in fields ranging from materials science to photonics, optoelectronics, and quantum technologies.

Typical setups for THz spectroscopy
To investigate the structure and optoelectronic properties of 2D materials using THz spectroscopy, various experimental setups have been established.The most common setup for THz spectroscopy includes a femtosecond laser for generating THz wave and a detection part based on electro-optic (EO) sampling.As we know the TES utilizes macroscopic photocurrents to generate detectable THz pulses that provide information about the dynamics of current flow.This technique has emerged as a crucial tool for investigating the behavior of photocarriers under femtosecond laser excitation in sophisticated material systems.Figure 2a shows a specific setup for TES 58 .A Ti: Sapphire oscillator generates ultrashort pulses that induce THz emission when incident on the sample at a 45°angle.Gold-coated parabolic mirrors collect, collimate, and focus the THz radiation onto a ZnTe (110) EO crystal.The electric field of the THz wave induces a change in the refractive index of the crystal, modulating the polarization of a gating beam passing through the crystal, and the gating beam is then measured using a differential detector consisting of two photodiodes.The THz electric field waveform is measured as a function of time by varying the time delay between the pump and gating pulses.It's important to note that a TES setup is a fundamentally broad concept, differing significantly from the specific setup mentioned previously.Depending on the systems to be studied, it can use various lasers with tuned wavelengths, and the measurement can be in either transmission or reflection geometry, and laser incident angle can be any angle, depending on the lattice structure and symmetry.The THz detector can use different nonlinear crystals depending on the THz intensity and bandwidth to be measured.TES has been employed in the investigation of semiconductors and extended to various other materials [59][60][61] .For instance, it was used to study aligned semiconducting carbon nanotubes, where the amplitude of THz emission correlated with the deposited pump energy and revealed the ultrafast dissociation dynamics of excitons 62 .The studies demonstrate the immense potential of TES in studying advanced material systems in the future.Moreover, TES serves as a valuable complement to THz time-domain spectroscopy (THz-TDS), through the spatial scanning of the pumplaser spot, it is possible to obtain localized information about the dynamics of photocarriers, achieving a spatial resolution limited only by the diameter of the pump-laser beam.The detailed exploration of TES, from its experimental setup to its diverse applications in studying material properties, underscores its significance in advancing our understanding of 2D materials and optoelectronics, paving the way for further innovations in material science.
In addition, THz-TDS is a widely used technique 63,64 , which utilizes ultrafast THz pulses for studying the interaction between THz radiation and matter.It involves measuring the temporal profile of THz pulses transmitted or reflected by materials to characterize their properties like complex refractive index.This non-destructive and contactless technique finds applications in diverse fields, including materials science, physics, biology, and telecommunications.A typical THz-TDS setup using photoconductive antennas or nonlinear crystals for both THz pulse generation and detection is shown in Fig. 2b 65 .A femtosecond infrared pulse is split into two, with one part generating the THz pulse and the other part used for detection via the EO sampling.This setup enables the characterization of materials in the THz range by measuring the transmitted THz waves for sample and reference.Specifically, photoconductive antennas (PCAs) are also widely used for detection, the process differs from EO sampling, though both methods aim to measure the THz electric field in the time domain.The THz-TDS allows phase-sensitive generation and detection of THz pulses, which directly provides the complex components of the material without the need for applying the Kramers-Kronig relations 66 , for example, magnetic susceptibility, optical conductivity, permittivity, etc., facilitating comprehensive investigations of different materials and phenomena.The THz frequency-domain spectroscopy (THz-FDS) system, utilizing continuous wave (CW) THz technology, offers precise frequency control and a higher dynamic range compared to the traditional THz-TDS, especially valuable for measurements requiring frequency sweeping and higher spectral resolution for sharp excitation probing.Figure 2c illustrates a typical THz-FDS setup, employing four-wave mixing, which achieves significantly better spectral resolution and dynamic range than the THz-TDS illustrated in Fig. 2b.In other aspects, by incorporating an optical pump path, the THz-TDS can be transformed into a so-called time-resolved THz spectroscopy (TRTS) technique, which enables the study of ultrafast carrier dynamics 67 .A specific TRTS system, depicted in Fig. 2d 68 , the pump beam is modified using a crystal to achieve the required wavelength and controlled fluence.THz probe pulses are created via optical rectification in a ZnTe crystal, and the THz radiation that passes through is detected using the EO sampling method in another ZnTe crystal.This setup allows for the investigation of the dynamics of photocarriers with high precision.By using TRTS, one can measure the changes in optical properties with sub-picosecond resolution, which provides insights into the relaxation dynamics of electronic states in the sample, facilitating the study of ultrafast quantum dynamics in materials 67,69,70 .
Based on TES, by integrating THz spectroscopy with other mature complementary characterization techniques, more comprehensive range of material characterization becomes achievable, enabling the exploration of local properties in a wider range of materials such as semiconducting nanostructures and 2D materials.For instance, the combination of THz spectroscopy with  58 .b, c Reprinted with permission from ref. 65 .d Reprinted with permission from ref. 68 .e Reprinted with permission from ref. 72 .f Reprinted with permission from ref. 78 scanning tunneling microscopy (STM) enables the attainment of high spatiotemporal resolution.The integration of THz-STM facilitates the study of ultrafast electron dynamics at sub-nanoscale spatial dimensions and sub-picosecond temporal resolutions [71][72][73] .This emerging technique has garnered significant attention due to its potential for uncovering intricate details and phenomena in low-dimensional materials.Figure 2e 74 .This work successfully visualized electron motion triggered by variations in the lowest unoccupied molecular orbital.This method will be crucial for evaluating local electronic structures and dynamics in 2D materials and semiconducting nanodevices.As the only technique which is capable of taking single-electron ultrafast movies at atomic resolution, THz-STM still has a lot of room for improvement, such as a higher signal-to-noise ratio 75 .Furthermore, the introduction of scattering-type scanning near-field optical microscopy (s-SNOM) has revolutionized subwavelength optics 76,77 .By coupling electromagnetic radiation to a subwavelength metal tip near a surface, s-SNOM allows the measurement of scattered radiation in the far field.Figure 2f shows the first s-SNOM measurement using a blue light (photon energy exceeds 3 eV) incident wave in the form of blue-pulse-induced THz emission 78 , high-energy pump photons enable strong THz emission from bulk crystalline silicon via two-photon excitation above the wide direct bandgap.Nanoscale spatial resolution is achieved in laser THz emission microscopy (LTEM) using femtosecond pulses at 410 nm and inducing THz emission with a sharp metallic atomic force microscope (AFM) tip.The setup involves generating near-infrared, blue light, and THz beams separately and detecting scattered or emitted THz pulses using EO sampling.In near-field measurements, the pump and probe spots are typically localized to the scattering probe's tip.However, to study non-local phenomena and obtain direct information about non-local effects, the use of two optical beams has been explored 75 .Pizzuto et al. conducted a non-local near-field pump-probe experiment on undoped bulk GaAs and GaAs nanowires.By photoexciting the sample at a laterally shifted location from the metallic probe tip using an ultrafast near-infrared pulse and probing with a broadband (0.2-2 THz) pulse scattered from the metal tip, researchers observed time-shifted conductivity responses consistent with carrier mobility simulations.This experiment confirmed the anisotropy of undoped GaAs nanowires and demonstrated the potential of the technique for studying nanoscale anisotropy, which is not feasible in the far field.The exploration of non-local phenomena through advanced THz spectroscopy techniques, as exemplified in the studies of bulk GaAs and GaAs nanowires, underscores the versatility and depth of THz analysis.The various THz spectroscopy and extended techniques mentioned facilitate comprehensive investigations into the structure and properties of materials.These advancements in THz spectroscopy open up new possibilities for studying ultrafast quantum effects, characterizing nanoscale structures, and exploring nonlocal effects, ultimately advancing our understanding of material physics and facilitating the development of advanced technologies.The coming sections will focus on the discussion of THz emission spectroscopy and review its application in studying 2D materials and their nanostructures.

THz emission from non-centrosymmetric 2D materials
Building on the foundation of THz emission spectroscopy's versatility as discussed earlier, we now delve into its application in the realm of 2D materials, particularly in the study of layered TMDs known for their exceptional optical properties.Layered TMDs, such as WS 2 , WSe 2 , MoS 2 and MoSe 2 , are the type of 2D materials widely known for their remarkable optical properties.Unlike gapless graphene, TMDs exhibit non-centrosymmetric lattice structure and possess tunable band gaps that can transition from indirect to direct band gaps at monolayer limit 31 .Figure 3a shows a typical experimental setup for THz emission from a layered 2D monolayer.The layered MoS 2 crystal generates efficient linearly polarized THz radiation (0.1-3.5 THz) under femtosecond laser excitation and the THz amplitude exhibits a linear dependence on pump fluence, indicating a second-order nonlinear process 79,80 .In principle, under femtosecond laser, charges in the materials are accelerated by the surface field in the surface-charge region.This acceleration leads to the generation of THz emission, which contains valuable information such as amplitude, polarization, and phase.Furthermore, THz surface emission spectroscopy enables non-contact observation of surface optical-physical processes with high temporal and frequency resolution.This technique has been successfully applied to investigate surface phenomena in various materials, especially in 2D materials.It is important to note that the THz emission mechanisms can differ between bulk materials and monolayer materials 81 .For example, the depletion field effect dominates in bulk WSe 2 82 , while the surface optical rectification dominates in bulk MoS 2 Commonly, optical rectification is one of the primary mechanisms in generating THz radiation using the second-order nonlinear effect.Taking MoSe 2 as an example (Fig. 3b), under an 800 nm pump excitation, the movement of positive and negative charge centers results in a change in the charge density, which leads to an alteration in polarization density 84 .Because the displacement of Mo and Se atoms occurs at the bond length level, and this process unfolds on a femtosecond timescale, the frequency of the excited electromagnetic wave falls within the THz range.Similarly, excitation of THz electromagnetic waves has been observed in other 2D materials, such as WSe 2 (Fig. 3c) 82 , indicating that optical rectification is mostly the dominant mechanism of THz emission in non-centrosymmetric monolayer materials.External conditions such as the polarization excitation in different directions of the pump can also impact THz excitation.For example, in WS 2 , when pump beams are in the s-polarization direction and p-polarization directions, they have an impact on the amplitude and phase of the excited THz radiation components, E x and E y .The amplitude of E x is greater than that of E y , indicating that the THz energy is mainly concentrated in the E x component, and the phase of E x and E y is roughly opposite (Fig. 3d) 83 .Therefore, circularly polarized excitation generates elliptically polarized THz waves (max.ellipticity ≈ 0.52) through circular optical rectification.In addition, polarized THz radiation can be easily controlled with the quarter-wave plate (QWP), which makes THz emission more controllable and advance the understanding of THz radiation properties and polarization control in 2D materials.The physical reasons for the different THz emission phenomena between monolayer and bulk materials are worth studying.).a, f Reprinted with permission from ref. 80 .b, e Adapted with permission from ref. 84 .Copyright 2024 American Chemical Society.Fig c, reprinted with permission from ref. 82 .d Reprinted with permission from ref. 83 .g, h Reprinted with permission from ref. 85 time-domain and frequency-domain spectra 84 .In the time-domain spectra, the THz amplitude of bulk materials is nearly four times that of monolayer materials due to the larger surface depletion field in bulk materials.In the frequency-domain spectra, the frequency range of radiation emitted by bulk materials is approximately 0-5 THz, while that of monolayer materials is 0-4 THz (Fig. 3e inset).For the bulk MoSe 2 , it has the point group structure of P 63 /mmc Space group and D 6h .This central symmetry will theoretically make it lose the second-order nonlinear effect.However, it is noted that the surface of the bulk material only has the periodicity in the XY direction, and the periodicity in the Z direction is broken, leading to the point group structure degenerating to D 3h , affecting the nonlinear coefficient leading to the generation of THz radiation.Although the surface of monolayer material and bulk material is similar, the width of the space charge region on the surface of bulk materials is much larger than that of monolayer materials.The presence of the space charge region induces the generation of transient photocurrent, and the first-order derivative of photocurrent on time determines the amplitude of THz radiation.Moreover, in bulk MoSe 2 , the influence of surface fieldinduced photocurrent on THz amplitude is greater than that of surface optical rectification on THz amplitude.Therefore, the amplitude of THz emission from bulk MoSe 2 is approximately four times greater than that of monolayer MoSe 2 .Additionally, 2D monolayer nanomaterials exhibit notable surface and interface effects, in-plane dipole effects, and strong spin-orbit coupling, leading to complex and intriguing nonlinear optical processes in THz emission compared to bulk materials.This analysis of THz emission mechanisms in both monolayer and bulk materials, particularly in the context of MoSe 2 and WS 2 , not only highlights the intricate interplay of physical processes governing THz radiation but also underscores the unique properties of 2D monolayer materials, setting a foundation for further exploration into their complex nonlinear optical behaviors and potential applications.
To explore the THz emission mechanism of different 2D materials, researchers have compared the THz emission graphite, MoS 2 , and InAs.In graphite, the presence of inversion symmetry in its layered structure prevents the occurrence of optical rectification.New mechanisms, such as the photon-drag mechanism, have been found to be responsible for its THz emission.For InAs, the carriers' effective mass and mobility are significantly different.The effective mass of holes is more than ten times that of electrons, and the hole mobility is less than 1% of the electron mobility, while the effective mass and mobility of electrons and holes in MoS 2 are in the same order of magnitude, which is common in 2D TMDs.This difference in InAs allows it to easily emit strong THz radiation by the photo-Dember effect, where the different mobility of electrons and holes leads to the destruction of surface symmetry, generating a Dember electric field.Under the modulation of the diffusion electric field and the Dember electric field, the excited electrons and holes form rapidly changing dipoles perpendicular to the surface and emit THz radiation.In addition, because the penetration depth of infrared light in MoS 2 is greater than that of the other two materials, there will be a period shift in the timedomain spectrum (Fig. 3f) 80 .In addition to the 2D materials mentioned above, researchers have also observed THz emission in quasi-2D hybrid organicinorganic perovskite (HOIPs) thin films via the photo-Dember effect (see Fig. 3g) 85 .The larger the laser photon energy, the larger the amplitude of the corresponding main peak of the THz wave.As shown in Fig. 3h, when the laser photon energy is greater than 2.43 eV and the time is 0.5 ps, the trend of these curves changes, the researchers illustrated two different polaron modes in the THz emission spectrum are existed, one dominated by inorganic sublattice vibration (P1) and the other involving the A-site cation (P2).This THz emission mode, which changes with the increase of excited photon energy, deserves further study.This exploration of THz emission mechanisms across a spectrum of 2D materials, including graphene, MoS 2 , InAs, and hybrid organic-inorganic perovskites, not only reveals the diverse underlying physical processes like the photo-Dember effect and photon-drag mechanism but also highlights the potential for further investigation into the complex interplay of material properties and THz emission characteristics, paving the way for new discoveries and applications in the field of nanotechnology.
Overall, the study of THz emission mechanisms in various monolayer (quasi-2D) materials provides insights into the complex processes involved in THz emission generation and helps us understand the unique properties of different materials in the THz regime.The explanation of these mechanisms allows researchers to unravel the 2D material physics underlying THz emission.In addition to investigating the THz emission mechanisms, there is ongoing research focused specifically on TMD-based THz emitters.These studies primarily aim to optimize and produce various types of controllable TMDs that can be utilized in THz applications.By engineering TMDs for specific THz applications, researchers can harness the unique properties of these materials and tailor them to meet the requirements of different technological advancements.We believe this ongoing research in TMDbased THz emitters holds great promise for the development of advanced THz devices and systems in the future.

THz emission originated from the photondrag effect
THz emission, mainly from the second-order nonlinear process, is commonly expected to happen in noncentrosymmetric crystals.However, in centrosymmetric materials, THz radiation can also be generated by introducing additional asymmetry, such as inversion symmetry broken by light-induced quasiparticles, surface effect, or photothermal effect, etc.Among these approaches, the photon-drag effect has been studied as an important mechanism which utilizes the incoming beam with a tilted incident angle.Such an effect has been used for the generation of THz emission in graphene.As sketched in Fig. 4a, the femtosecond optical pump pulse illuminating at oblique incidence induces nonzero shift current dipole arising from an asymmetric distribution of nonthermal electron and hole population even in centrosymmetric multilayer graphene, giving rise to THz emission covering 0.1-4 THz 86 .This nonzero shift current dipole arises from a nonvertical transition between the valence and conduction bands due to the transferring of the finite inplane photon momentum to the electron-hole pairs (Fig. 4b) 87 .It is worth noting that the polarity of the signals from graphene reverses for opposite incident angle (Fig. 4c) 88 , corresponding to an opposite shift current dipole caused by the opposite in-plane photon THz peal-valley value (arb.units)  86 .b Adapted with permission from ref. 87 .c Reprinted with permission from ref. 88 .e Reprinted with permission from ref. 89 .f Reprinted with permission from ref. 90 .g Reprinted with permission from ref. 91 .h Reprinted with permission from ref. 92 .i Reprinted with permission from ref. 93 .j Reprinted with permission from ref. 94 .k Reprinted with permission from ref. 95 momentum.This relationship is further confirmed by varying the incident angle and comparing the peak-topeak amplitude of the emitted THz electric field normalized by (1 − r) (r is the amplitude Fresnel coefficient of graphene) (Fig. 4d) 86 .Such mechanism based on photon-drag is similarly validated in another work with the surface photogalvanic effect excluded, in which the polarization of the laser beam can modulate the THz emission intensity 89 .Concretely, the excitation using a p-polarized beam from the substrate side reverses the polarity of the generated THz wave (Fig. 4e).
For improved efficiency and extra tunability of the photon-drag effect, THz emission from vertically grown graphene (VGG) is achieved.Figure 4f compares the emitted THz signal from VGG with monolayer graphene, the enhanced signal is attributed to the light-matter interaction enhancement and the more intensive photoexcited carriers in the VGG structure 90 .Additionally, the enhanced efficiency also enables prominent detection of circular photon-drag effect (CPDE).Helicity-dependent THz emission is observed in VGG by tuning the QWP polarization angle 91 , which is confirmed as the result of the combined contribution of CPDE and linear photondrag effect (LPDE) with different weights As exhibited in Fig. 4g, there is a regular cycle of change in the peak-valley values of THz p components in the form of cos 4φ (φ is QWP polarization angle) due to the dominance of CPDE, while the amplitudes of THz s-components exhibits irregular variation curves as the weight of CPDE and LPDE changes with φ.
Apart from VGG, THz emission with modulated amplitude and phase is also demonstrated in layered selenides like Bi 2 Se 3 and SnSe 2 .As an illustration, the polarized light excitation in the Bi 2 Se 3 film could induce photon-drag effect (Fig. 4h) 92 , in which an apparent periodicity of the azimuth angles is observed (shown in the inset), attributed to the intrinsic symmetry of the Bi 2 Se 3 crystal.Further research on SnSe 2 demonstrated that applying a constant electric field is an effective way to enhance the intensity of photon-drag induced THz radiation to 123% (Fig. 4i) 93 .
Recently, PtSe 2 attracted attention as an emerging 2D platform with a centrosymmetric crystal structure, in which THz emission could be observed when the incident angle (θ i ) of the pump pulses is not zero (Fig. 4j) 94 .The emission efficiency through photon-drag mechanism is reported up to 5.2 × 10 10 (defined as E/Fd, where E is the maximum THz electric field, F is the pump fluence, and d is the penetration depth for 800-nm light or the thickness of the sample), as a result of its particle-hole asymmetric band structure.For more specific applications, such THz emission from PtSe 2 with polarization dependence was demonstrated to realize ultrafast ternary encoding as a new detection model (Fig. 4k).The positive or negative polarities of THz pulse signals were generated by leftcircularly polarized or right-circularly polarized laser excitations, which was coded as 1 or −1, respectively.In contrast, when using a linearly polarized laser to excite, no THz pulse signals were detected which was coded as zero 95 .Such THz coding function based on photon-drag effect spreads the non-contact modulation capability and the application prospect in THz emission applications.

THz emission from heterostructures
2D materials have the advantage of easy integration with traditional optoelectronic materials, systems, and devices, leading to novel phenomena and applications, especially in the THz regime.Integration is achieved through direct growth or thin-film transfer methods.The heterostructures can be usually distinguished into two types: one for vdW heterostructures are composed of different 2D materials 96,97 and the other is heterogeneous integration with traditional materials and devices 98 .
As one of the typical 2D materials, Graphene has long been used to construct heterostructures to emit THz.It is commonly combined with SiO 2 and Si bases for the onchip industry.In a study by Yao et al., a noninvasive and efficient method was applied to investigate electric-fieldinduced optical rectification that indicates the correlation between built-in interface potential and THz generation.Without applying a gate voltage, Gr/SiO 2 /Si generates stronger THz than Si (Fig. 5a) 54 , which is mainly from the photocarriers with the picosecond response in the interface layer.When applying a negative gate voltage, THz generation is enhanced, where THz amplitude is three times higher when V g = −2V compared with V g = 2 V (Fig. 5b).This indicates that modifying the interface depletion layer can tune THz generation with gate voltage applied.Graphene/Si Schottky junction (GSSJ) is the efficient THz interfacial emission structure under the gate voltage.
Similar to monolayer materials, the polarization of heterojunction THz emission has also been studied.The orthogonal components of the induced THz field include p-and s-polarized THz emission, with the amplitude of s-polarized emission being much lower than that of p-polarized THz radiation (Fig. 5c) 99 .The investigation of THz emission in heterostructures related to TMDs has garnered significant attention.Compared with WSe 2 only, WSe 2 /Si THz emission amplitude is ~10 times larger under a 1.65 eV laser pump (Fig. 5d) 98 .Additionally, THz modulators based on other TMDs-Si heterostructures, such as MoS 2 -Si, have also been fabricated and investigated.Since MoS 2 sheets are nearly transparent for THz waves, THz emission amplitudes of MoS 2 and Si metasurfaces are nearly the same when no laser is pumping.Compared to the reference signal in air, the pulse broadening and deformation are attributed to the resonances within the metasurface structure.Under laser pumping of 1.75 W/cm 2 , the amplitudes of THz emission from both metasurfaces decline, and the amplitude of THz emission from MoS 2 decreases even more.This phenomenon can be attributed to the enhanced modulation depth of the Si/MoS 2 configuration, leading to more photocarrier generation compared to a purely Si surface (Fig. 5e) 55 .
The heterostructures combined with TMDs and other 2D materials have also been widely studied in the field of THz emission.Krihna et al. applied hBN/WS 2 heterostructure to investigate THz photoconductivity and photocarrier dynamics with interfacial dipole formation at the interface between insulator and semiconductor 96 .Other vdW heterostructures such as ReSe 2 /MoS 2 heterostructure have also been investigated using THz emission (Fig. 5f inset).As shown in Fig. 5f, THz radiation for both ReSe 2 /MoS 2 and MoS 2 /ReSe 2 direction last for 300 fs 100 .In this study, various stacking orders and measurement techniques were employed to investigate photocarrier dynamics.Reversal of the emitted THz pulse was observed when the stacking order was interchanged, indicating that THz emissions arise from the plane current, and the frequency spectrum of the waveforms exhibited a peak at 1.0 THz with a bandwidth of 2.5 THz.A similar phenomenon has also been observed in  54 .c Reprinted with permission from ref. 99 .d Reprinted with permission from ref. 98 .e Reprinted with permission from ref. 55 .f Adapted with permission from ref. 100 , Copyright 2024 American Chemical Society.g Reprinted with permission from ref. 56 .h Reprinted with permission from ref. 101 .i Reprinted with permission from ref. 103 .j Reprinted with permission from ref. 104 .k Adapted with permission from ref. 97 .l Reprinted with permission from ref. 105 WS 2 /MoS 2 heterostructure that emits much stronger THz emission than pure WS 2 or MoS 2 , indicating the interfacial charge transfer dominates the THz emission process 56 (Fig. 5g).Perovskite heterostructures, going beyond the realm of common heterostructures, have garnered significant interest in scientific research.These perovskite devices are particularly notable for their high modulation efficiency, rapid response times, cost-effectiveness, ease of integration, and stability.Specifically, the amplitude spectrum of the CsPbBr 3 perovskite quantum dots combined with Si under various laser intensity irradiances is illustrated in Fig. 5h 101 .Additionally, the spectral transmission of the CsPbBr 3 perovskite quantum dot heterostructure was meticulously measured using a Backward Wave Oscillator (BWO).Although it is not based on TES but rather TRTS method, there are still abundant physical phenomena worth exploring.These measurements were taken both with and without the application of an external 450 nm pump laser excitation, across a frequency range from 0.23 to 0.35 THz.To delve deeper into the modulation mechanism, researchers analyzed the optical constants of the CsPbBr 3 perovskite quantum dots/Si.This analysis was conducted through time-frequency spectral measurements using the BWO under a range of laser irradiances.The findings revealed that the CsPbBr 3 perovskite quantum dot heterostructure boasts a modulation speed of 2.5 MHz and a modulation depth of 45.5%.These impressive performance metrics strongly suggest that this heterostructure could be an excellent candidate for future applications in practical THz wave communication systems.
In addition to 2D materials, heterostructures composed of traditional materials also have great potential in THz emission.Kampfrath et al. demonstrated that laser pulses irradiating heterostructures composed of a ferromagnetic iron layer and a non-magnetic metal (ruthenium or gold) layer drive ultrafast spin currents and generate THz emission.This discovery indicates the great potential of spintronic THz emitters (STE) as strong THz sources 102 .Based on these findings, other heterostructures are also being utilized for THz emission.For example, Co/Bi(n)/Bi 2 Te 3 multi-material stacking heterostructure is found to emit strong THz radiation.Adding (or reducing) a certain layer of material can cause a sharp change in the efficiency of heterojunction THz emission.Compared to the negligible THz signal from Bi 2 Te 3 , Bi/Bi 2 Te 3 , or Co/Bi films, strong THz emission is found from Co/Bi 2 Te 3 , Co/Bi/Bi 2 Te 3. due to the spin-to-charge conversion process (Fig. 5i) 103 .Significantly, the amplitude of the THz signal emitted from the Co/Bi/Bi 2 Te 3 heterostructure can reach up to 198% of that emitted by devices lacking an intercalated Bi film.
Notably, shift current is another strategy to emit THz radiation in the heterostructures.In the GaAsSb 0.13 /In 0.141 Ga 0.859 As heterostructure, lower pump-pulse duration leads to an increase of the peak field of THz transients from 0.51 to 1.1 kV cm −1 .On the other hand, the width of the field crest Δτ decreases with the pump-pulse duration from 70 fs to 20 fs (Fig. 5j) 104 .In addition to the above-mentioned combination of Graphene and traditional materials to build heterostructures, Graphene also builds heterostructures with other 2D materials.For example, in 2020, graphenelayer/black phosphorus (GL/BP) and GL/MoS 2 heterostructures are theoretically expected to be THz emitting diodes (THz-EDs) with lateral hole and vertical electron injection.The increasing injection current density may lead to stronger output power in both heterostructures (Fig. 5k) 97 .The increase is correlated with doping level (μ a ), and undoped GL is shown to be significantly boosted by injection current.It is noteworthy to clarify that, in this research, the term "injection current" pertains to an electrical means of THz emission, distinct from the optical counterpart.Although further experimental verification needs to be completed, this work has provided a feasible method for fabricating the more advanced THz emission devices based on 2D materials.
Researchers are increasingly focusing on complex heterostructures, particularly those incorporating a variety of materials.A notable example is as shown in Fig. 5l 105 , the CH 3 NH 3 PbI 3 /graphene heterostructure, where under an 800 nm, 288 mW pump beam, the THz peak signal initially rises then falls to a negative value, indicative of photoenhanced transparency similar to monolayer graphene's hot-carrier effect.This effect increases graphene's chemical potential and temperature, while carrier−carrier scattering reduces electrical conductivity.In the CH 3 NH 3 PbI 3 /Ag wire grid (Ag WG)/graphene sample, the Ag wire grid boosts pump beam absorption, enhancing photogenerated carrier efficiency and transient THz transmittance.This sample's THz signal recovers quicker than single-layer graphene's, due to faster hot-carrier cooling and photoexcited carrier accumulation in the CH 3 NH 3 PbI 3 layer.Furthermore, the sample exhibits a delayed negative THz signal, attributed to the long carrier migration distance in the CH 3 NH 3 PbI 3 layer and the transfer of hot carriers from graphene to perovskite, thereby enhancing photoconductivity.The total THz transmission reflects the interplay between graphene and CH 3 NH 3 PbI 3 materials.The Ag wire grid aids in electron transfer in CH 3 NH 3 PbI 3 , promoting effective electron-hole separation and efficient charge collection on graphene.This innovation opens new ways for THz wave modulation and visible light detection applications.
In summary, the investigations of THz emission from heterostructures have revealed their remarkable potential for generating and controlling THz radiation.Whether traditional material heterostructures or vdW heterostructures, these materials have demonstrated efficient THz emission through optical rectification, interfacial dipole formation, shift current, and spin-to-charge conversion which will be discussed in detail in the following section.The manipulation of gate voltages, modulation depths, and shift currents allows for tunability and enhancement of THz generation.Moreover, utilizing heterostructures opens up possibilities for practical applications and offers a versatile platform for exploring and harnessing the potential of THz radiation in various fields, paving the way for advancements in THz technology.

High-order harmonic THz generation
Distinguished from the typical THz emission induced by optical rectification and photon-drag effect, the THz high-order harmonic generation (HHG) reflects the nonlinear process beyond second-order in materials including graphene, TMDs, Dirac semimetals, and topological insulators.For instance, an incident quasi-monochromatic, linearly polarized pump wave of 0.3 THz will cause harmonic emissions up to the seventh-order from the monolayer graphene at room temperature and under ambient conditions (Fig. 6a) 106 .In principle, the multiplication of photon energy is driven by the nonlinear interaction between light and matter.Free background Dirac electrons in graphene feature nonlinear intraband THz conductivity originated from the ability of internal thermalization and plays a role as an energy reservoir that promotes the light-to-matter energy transfer process.The nonlinear conductivity consequently results in modulated THz absorption and current in response to the pumping wave, leading to THz harmonic re-emission at multiplicated frequencies.On account of the centrosymmetric nature of graphene, only odd THz harmonics are generated.Moreover, parameter-free calculations using a thermodynamic model verify the high-harmonic generation efficiency and illustrate the potential for up to 13th THz harmonics (Fig. 6b) 106 .Particularly, the key to highefficiency HHG emission is introducing and regulating the background carriers.By applying a gated voltage to graphene, the Dirac point and the background carrier density can be adjusted to the appropriate value, thus achieving the best HHG emission efficiency by striking a balance between enhancing the power absorption of the THz drive field and avoiding the nonlinear thermodynamic decline (Fig. 6c) 107 .Through this method, an increase of about two orders of magnitude in transmission efficiency has been achieved (Fig. 6d) 107 .
In addition to graphene, THz high-order harmonic generation is also observed from other Dirac semimetals such as Cd 3 As 2 .With the nonlinear kinetics of electron distribution and linear energy-momentum dispersion, electrons in the Dirac cone are accelerated and scattered by THz electric fields, accumulating energy into the electron subsystem, leading to a stretched and shifted distribution along the field, resulting in efficient HHG (Fig. 6e) 108 .Specifically, the HHG from Dirac semimetals exhibits non-perturbative fluence dependence and is predicted to be highly sensitive to the scattering rate.Nevertheless, the HHG characteristics of Cd 3 As 2 are reported to be related to the supercritical or subcritical regime of nonlinear optics, corresponding to whether all electrons can be driven beyond the Dirac point 109 .In the supercritical regime, up to the 31st order HHG can be generated using a modestly powered driving laser (Fig. 6f).In contrast, the subcritical regime leads to zero intraband emission of 5th and higher order HHG.More recently, the regulation of third-harmonic yield and polarization state in Cd 3 As 2 through driving-pulse ellipticity has been achieved (Fig. 6g) 110 .The increased ellipticity results in enhanced perpendicular component and decreased parallel component of the emitted electric field through nonlinear intraband kinetic processes, indicating the THz optical signal processing potential of Dirac semimetals.
Besides, as a type of Dirac material, topological insulators (TI) such as Bi 2 Te 3 exhibit interesting THz HHG characteristics.Innovatively, tunable non-integer HHG is demonstrated in Bi 2 Te 3 film 111 .With the ballistic acceleration of Dirac currents in the topologically protected surface state, HHG dominated by the bulk or the surface state can be selected by tuning the driving frequency (Fig. 6h).When the carrier-envelope phase of the driving field changes, the special dynamics enables a continuous shift of HHG to arbitrary non-integer multiples of the driving frequency.Furthermore, topological insulators demonstrate intriguing feature of enhanced THz HHG efficiency.In topological insulator Bi 2 Se 3 , third-order THz harmonics at milliwatt level with a field conversion efficiency of ~8% is achieved, which is nearly two orders of magnitude better than monolayer graphene at a high incident power (Fig. 6i) 112 .The enhanced emission is attributed to the ultrafast dissipation of electronic heat through the surface states of the bulk and thus avoiding heat accumulation in both electronic and phononic systems.Moreover, using a metal grating combined with TI can increase the harmonic power from low to intermediate incident powers by surface-selective facilitating.Apart from TI, THz HHG can also be generated in 2D semiconductors such as TMDs.In MoS 2 , up to the 13th-order harmonic can be observed when excited by incident light of 0.3 eV (72.55 THz) 113 .Due to the non-centrosymmetric nature of its lattice, even-order harmonics also appear (Fig. 6j).Compared to bulk MoS 2 materials, the enhanced odd HHG efficiency is considered relative to the reduced dielectric screening and strong electron-hole Coulomb interaction in the monolayer.
Clarification of the nonperturbative electron dynamics under low-frequency and high-intensity fields in these Dirac systems is of great significance.Recently, an analytical approach based on the (Jeffreys-)Wentzel-Kramers-Brillouin (WKB) approximation has been proposed to analyze HHG in the THz regime for realistic massive Dirac materials, which provides HHG waveforms that match well with the numerical results (Fig. 6k) 114 .THz HHG in the material system with Dirac electronic properties provides additional emission frequencies and offers valuable insight on the intraband dynamics of Dirac electrons.However, the frequencies and relative intensities of the harmonics are usually limited to multiples of the excitation frequency and magnitudes of the incident power several orders below, which requires further exploration to expand its availability for THz source applications.

THz emission arising from spin-to-charge conversion
Different from the mechanism that requires the symmetry-breaking either by the non-centrosymmetric lattice structure or by the asymmetrical geometry of the photoexcitation, a novel mechanism known as spin-tocharge conversion has emerged as an alternative pathway for THz emission in heterostructures comprising a solitary ferromagnetic (FM) layer alongside a non-magnetic (NM) layer.In such heterostructures, the spin-orbit coupling at the FM/NM interface plays a crucial role in the spin-to-charge conversion process.In brief, while an ultrafast laser pulse is incident on the heterostructure, it generates a non-equilibrium distribution of spinpolarized charge carriers in the FM layer, wherein the spins of these carriers are coupled to their momentum due to the spin-orbit coupling at the FM/NM interface.This coupling leads to the generation of a pure spin current in the FM layer, which propagates across the FM layer until it reaches the FM/NM interface.At this interface, the spin current can undergo a conversion into a charge current, leading to the emission of THz radiation.In essence, this spin-to-charge conversion process occurs due to the transfer of angular momentum between the spin and the charge degrees of freedom at the interface.
Such spintronic THz emission is observed from metallic FM/NM heterostructures in 2013.When Fe/Au k HHG spectrum for the oscillating field with analytical WKB theory.a, b Reprinted with permission from ref. 106 .c, d Reprinted with permission from ref. 107 .e Reprinted with permission from ref. 108 .f Reprinted with permission from ref. 109 .g Reprinted with permission from ref. 110 .h Reprinted with permission from ref. 111 .i Reprinted with permission from ref. 112 .j Reprinted with permission from ref. 113 .k Reprinted with permission from ref. 114 or Fe/Ru heterostructures are incident with a femtosecond laser pulse at 1.55 eV, the absorption in the Fe layer made the majority-spin electrons in iron with mainly splike character excited, the non-equilibrium electron distribution then leads to immediate spin current (Fig. 7a), resulting in THz emission through the inverse spin Hall effect 102 .The THz emission polarizes perpendicularly with respect to the magnetization, covering the bandwidth from 0.3 ~20 THz, which is tunable by adjusting the non-magnetic layer material (Fig. 7b).Especially, excited by circularly polarized optical pulses, similar FM/ NM heterostructures composed of NM material with strong spin-orbit interaction exhibit helicity-dependent THz emission characteristics (Fig. 7c).For example, Co/Pt heterostructure can generate laser-induced THz radiation parallel to the magnetization M (E y ) 115 , which is distinguished from the pump polarization independent spintronic emission induced by the inverse spin Hall effect (E x ).The maximum amplitude of E y is typically one to two orders of magnitude smaller than E x , and hardly change with the thickness of the metal layer, indicating a spin-orbit interaction at the contact interface.Notably, the sign of the emitted THz radiation will change if either the helicity of the circularly polarized light reverses or the heterostructure is rotated around the magnetization over 180°(Fig..c Reprinted with permission from ref. 115 .d Reprinted with permission from ref. 116 .e, f Reprinted with permission from ref. 117 .g Reprinted with permission from ref. 118 .h Reprinted with permission from ref. 119 .i Reprinted with permission from ref. 120 .j Reprinted with permission from ref. 121 .k Reprinted with permission from ref. 122 .l Reprinted with permission from ref. 123 The output of the STE can be by designing the layer materials and multilayer film structures.For the FM/NM metal heterostructure, material of the NM layer significantly alters the THz amplitude, as shown in Fig. 7d.Among the various NM materials studied, the output intensity is not strictly correlated with the material's spin Hall conductivity 116 .Interestingly, W with a half-filled delectron shell results in a relatively large amplitude in an opposite direction.Exploiting this, placing two NM thin films with opposite spin Hall angles (W and Pt for instance) on either side of the FM layer enables the forward and backward spin currents to eventually convert into in-phase radiation, significantly improving the efficiency of pump energy utilization.As for the FM layer, the influence on the output amplitude is relatively small when different materials are used, among which the CoFeB film has been identified as a more optimal solution.Ultimately, the optimized W/CoFeB/Pt three-layer spin emitter can generate pulses covering the 1-30 THz range without a gap, with amplitudes surpassing standard emitters such as ZnTe, GaP, and biased photoconductive switches.
Apart from metal heterostructures, the challenge of giant spin injection into semiconductors is resolved in a Co/2D semiconductor MoS 2 heterostructure, leading to THz emission induced by spin-to-charge conversion 117 .As the out-of-equilibrium carrier distribution in the ferromagnetic layer diffuses and relaxes in energy states, almost only the majority-spin electrons remain at high energies, which are then filtered by the band gap of MoS 2 and result in a large spin injection effect (Fig. 7e).Throughout this process, the strong spin-orbital interaction of MoS 2 enables efficient THz emission induced by spin-to-charge conversion (Fig. 7f).Particularly, full parameter-space characterization shows that the THz amplitude has a spin-current-based origin, and is independent of the orientation-related nonlinear optical response.The peak THz intensity is regulated by pump photon energy, indicating a non-thermalized process in mechanism.Additionally, spin-to-charge conversion is realized in Bi 2 Se 3 /Co heterostructure, which is attributed to the important role of spin-momentum locked surface states in the topological insulator 118 .THz emission of the heterostructure is very apparent compared with the emission of either Bi 2 Se 3 or Co film (Fig. 7g), indicating a high efficiency of spin-injection and spin-to-charge conversion.The special properties of 2D materials can also be exploited in the ferromagnetic layer.2D vdW ferromagnetic Fe 3 GeTe 2 (FGT) has been used to construct FM/TI heterostructures 119 .FGT integrated with Bi 2 Te 3 exhibits spin-to-charge conversion-induced THz emission through the inverse Edelstein effect mechanism (Fig. 7h), regulated by the direction of the magnetic field, the pumping incidence, and the thickness of the two films.It is noteworthy that the interfacial exchange coupling in the FGT/TI heterostructure is observed to stabilize the ferromagnetic order of 2D FGT, providing insights into lowdimensional characteristics of layered FM materials.
For device applications in THz communication, modulation of amplitude, phase, polarization, etc. is significant for optimized transmission efficiency.Inspired by the linear polarization of spintronic THz emission from FM/ NM structure, a nonuniform external magnetic field is designed to realize THz wave generation and tuning of the chirality (Fig. 7i) 120 .The azimuthal angle can be manipulated by rotating the sample and the magnets.Additionally, the ellipticity can be tuned by changing the phase difference or the emitting areas, and also depends on the frequency.However, strict tuning of the phase difference for circularly polarized emission is still challenging, requiring further exploration in adjusting the material thickness and conductivity of the FM/NM nanofilms.
Coherent control of THz radiation is demonstrated in 2D hybrid metal halides (2D-HMH)/NiFe heterostructures 121 .As an emerging type of synthetic semiconductor, 2D-MHM features spin-orbit coupling due to the periodically layered inorganic PbI 6 framework structure, revealing the intriguing potential for STE.Through an inverse Rashba-Edelstein effect mechanism extended with additional in-plane momentum shift of the Rashba bands, the phase and intensity are manipulated by tuning the applied magnetic field and linear polarization of the femtosecond laser pump pulses (Fig. 7j).This dependence is not observed in Pt/NiFe or 3D-HMH/NiFe devices, indicating a low-dimensionality related origin of the 2D-HMH layer under below-gap femtosecond laser excitation.Moreover, by integrating Pt/CoFe heterostructure with a single crystal piezoelectric Pb(Mg 1/3 Nb 2/3 ) O 3 -PbTiO 3 (PMN-PT) layer, an artificial magnetoelectric coupling is endowed through strain transfer across the ferroelectric (FE)/FM interface, thus realizing an effective electric-field modulation of 270% of the THz amplitude (Fig. 7k) 122 .
In addition to exploiting material properties, surfacepatterned structural design is used for spintronic THz modulation and coded emission.Combining interspaced Co/Pt and Co/W heterostructures with an optical mask enables selective excitation of micro emitters with opposite spin Hall angles, resulting in coded THz phase components (Fig. 7l) 123 .Such spatial design provides an additional operational dimension, demonstrating the direction of more integrable and multifunctional devices.
Studying THz emission from heterostructures with FM and NM layers allows researchers to explore spin dynamics, spin transport, and spin-related phenomena.Meanwhile, the inclusion of extra spin degrees of freedom and the utilization of diverse FM/NM material systems broaden the possibilities for generating and modulating THz emissions, offering potential applications in spin-THz devices.

Extended THz sources and their applications
Beyond THz generation induced by femtosecond optical excitation, the utilization of various mechanisms and structures is also being thoroughly explored.For instance, the Smith-Purcell effect arises when charge carriers move on periodically patterned surfaces, leading to THz emission.However, traditional SP emitters are constrained by electron beam acceleration and alignment.Recently, an SP THz emitter utilizing charge carriers in 2D materials has been designed 124 .Numerical simulations indicate that coherent THz radiation can be expected by integrating graphene and silicon grating and exciting the hot electrons within (Fig. 8a).The performance of this emitter holds the potential for further modulation through experimental transformations of expanded 2D materials and the design of specific metasurface structures.It represents a promising THz source with potential for high efficiency, near-field excitation, and on-chip integration adaptability.
In comparison to THz time-domain systems, QCLs exhibit more availability as miniaturized direct laser sources, which provides a compact solution with chipscale dimensions, high power emission, high spectral purity, and wide bandwidth.However, engineering intracavity semiconductor multilayers with the necessary gain and absorption properties at THz frequencies presents significant challenges due to the extremely low photon energies involved.Riccardi et al. 125 demonstrated the generation of 4 ps pulses in a semiconductor heterostructure THz laser using passive mode-locking.A distributed graphene saturable absorber (DGSA) was produced by transferring a seven-layer CVD-grown multilayer graphene (MLG) film along the entire cavity of double-metal QCL (Fig. 8b).Commercial femtosecond systems in the visible/infrared spectrum rely on passive pulse generation due to their ability to generate the shortest and most powerful pulses.Consequently, this research offers promising prospects for spectroscopy applications.
Due to its unique frequency characteristics, THz light sources exhibit distinct advantages in fields such as ultrafast imaging, biosensing, and high-bandwidth communication.Leveraging the penetrative nature of THz radiation, a single-shot ultrafast imaging system has been realized 126 .The initial probing beam is divided into multiple sub-pulses differentiated in the frequency domain with varying delay times, thus achieving timeresolved image acquisition (Fig. 8c).Through the composite detection, decoding, and reconstruction in both  124 .b Reprinted with permission from ref. 125 .c Adapted with permission from ref. 126 .d Reprinted with permission from ref. 127 .e Reprinted with permission from ref. 128 .f Reprinted with permission from ref. 129 .g Reprinted with permission from ref. 130 .h Reprinted with permission from ref. 15 time and spatial frequency domains, this system possesses the capability capture transient events within nontransparent media, achieving a time resolution at the subpicosecond scale.This method extends the limitations of ultrafast photography based on traditional wavelength light source, providing a powerful tool for exploring fundamental dynamic events in various domains including physics, chemistry, biology, materials, etc.Moreover, implementing label-free biosensors in the THz spectra has enabled many potential applications, such as bio-material detection in point-of-care (PoC) diagnostics.Bai et al. 127 conducted a study in which they exhibited a monolithic THz emission biosensor integrating asymmetric double-split ring resonator metamaterials with a ferromagnetic heterojunction STE (Fig. 8d).This THz emission biosensor has enormous potential for highresolution near-field biosensing applications, especially for detecting trace biological samples.
The advancement of THz-frequency waveguides is necessary for developing integrated, cost-effective, and efficient solutions for high-speed devices.Yang et al. 128 demonstrated the robust transport of THz topological valleys through multiple sharp bends on an all-silicon chip to surmount these obstacles using full-wave simulations and visualized the intensity distribution of the magnetic field near the twisted domain wall (Fig. 8e).The THz topological photonic circuit is immune to sharp bends, and is able to achieve error-free transmission of uncompressed 4 K high-definition video, demonstrating the immense potential of THz topological photonics in applications such as topological splitters, robust delay lines, compact interferometers, and directional antennas.
To address integrated application scenarios, THz emission or detection can be described in the logic way of electronic and photonic way, which serves as the basis for THz technology.Figure 8f 129 demonstrates that TSMC has accomplished a remarkable feat by effectively fabricating a circuit using standard 65 nm CMOS technology with a die size of 1.55 × 1.80 mm 2 .Several crucial techniques, including incorporating band-selection inductor switches, simultaneous bulk voltage tuning across oscillators and dividers, a mutual injection locking mechanism, and a dual port injection divider, have been implemented to achieve a locking range of 21 GHz.
Significant advances have been made in developing THz signal generation technologies over the recent years (Fig. 8g).In the field of electronics, III-V-based semiconductor technologies have demonstrated the capacity to generate 100 μ W to mW range of power at room temperature 130 , representing a two-to five-fold increase, especially at frequencies above 1 THz.Similarly, photonic approaches to THz generation have made significant strides.Photomixer-based methods have achieved continuous-wave power levels in the tens of microwatt range.In contrast, difference-frequency or parametric generation utilizing nonlinear materials has achieved power levels in the 1-2 THz range, reaching the milliwatt range.
In the quest for extending the scale of THz phenomena investigation beyond the diffraction limit, significant attention is directed towards the synergistic integration of THz radiation with scanning near-field optical microscopy (SNOM).This technique is commonly categorized as d-SNOM assisted by subwavelength probe, and s-SNOM enhanced by a sharp tip.At different spatial resolutions and scales, the application space of near-field THz imaging includes bio-med imaging, THz devices, physics and materials, etc. (Fig. 8h).Stem from the long wavelength nature of THz radiation, the confinement of THz radiation to sub-wavelength scales and accomplishing nonlinear near-field spectroscopy in the THz range are the main challenges unresolved 15 .Breakthroughs in these limitations will propel advancements in THz research, particularly in the areas of cryogenic temperatures, ultrafast timescales, and nanoscale lengths.

Summary and future perspectives
To encapsulate, THz emission spectroscopy has proven to be a powerful and versatile tool for investigating the electronic and optical properties of 2D (or quasi-2D) materials and their interfaces.This review has explained the diverse nonlinear processes underlying THz emission, including optical rectification, photon-drag effect, interfacial processes, high-order THz harmonics generation, and spin-to-charge conversion, showcasing the rich diversity of THz emission spectroscopy in unraveling the complex nature of these materials.These findings contribute to the understanding of novel 2D (or quasi-2D) materials and pave the way for advancements in fields ranging from materials science to photonics, optoelectronics, and quantum technologies.From the perspective of technological application, these 2D (or quasi-2D) materials have shown huge potential to be nextgeneration laser-driven table-top THz emitters due to their atomically thin structure and easy integration into silicon wafer to make compact devices.Although these different 2D material and their structures can emit THz radiation under fs laser illumination, the real application using them for THz spectroscopy demands high efficiency and broad bandwidth to extract the THz response over a large range of spectrum.Of all the materials and structures mentioned earlier, STE stand out with significant advantages in terms of both high efficiency and a wide spectral range.Several strategies can be adopted to design the STE.First of all, the bilayer structure FM/NM structure can be replaced by a trilayer NM 1 /FM/NM 2 structure 116 , where the NM 1 and NM 2 have opposite spin-hall angles.Secondly, multilayer structure (NM 1 /FM/NM 2 ) n can be made to enhance the THz emission.Thirdly, most of the femtosecond laser pulse are wasted due to the transmission and at interfaces, whereas photonic crystal can be used to confine the laser electric field and strengthen the interaction between the laser and the metallic layer 131 .For example, constructing a simple Fabry-Pérot cavity is proved effective.The STE is sandwiched between a MgO layer and dielectric mirror [SiO 2 (165 nm)|TiO 2 (94 nm)] 5 , which confines the optical filed inside the cavity and minimizes pump transmission into Si and reflection into air 131 .Fourthly, special design such as THz anti-reflection film can be used to further enhance of the THz transmission.These strategies can be utilized either individually or in combination to attain maximum THz emission in the far field.Building upon this, explorations including the introduction of diverse 2D materials with strong spin-orbit interactions 132 , along with the complex structural control of the emitted THz waves, represents a frontier exploration in this field 133 .
However, there are still remaining challenges that need to be addressed to fully unlock the potential of THz spectroscopy of emerging 2D materials and heterostructures.This summary investigates the current problems and future development trends of THz spectroscopy, with a focus on achieving Spatial resolution, broadband spectral coverage, high signal-to-noise ratio, high efficiency and integrability.First of all, one primary challenge THz spectroscopy is spatial resolution, as many emerging materials, particularly 2D materials, often exhibit nanoscale features that are significantly smaller than the diffraction-limited spot size of the farfield THz probe (~300 µm).As a result, far-field THz spectroscopy is limited to spatially-averaged measurements, providing information across inhomogeneous regions or ensembles.Despite the availability of nearfield techniques discussed before, it remains difficult to achieve nanometer-scale resolution and probe individual nanoparticles or the internal structure of materials.Secondly, it is important to achieve efficient broadband coverage for THz generation and detection, as traditional THz sources and detectors often have limited frequency ranges, the necessity for wide-ranging measurements spanning from 0.1 to 30 THz is evident, as it covers the extensive frequency range encompassing spectral signatures in emerging materials.Additionally, the spectral characteristics may vary with experimental conditions and material geometry as discussed before.The THz conductivity spectra of emerging materials are frequently broadened due to high charge carrier scattering rates, further emphasizing the requirement for broadband sensitivity.To address this challenge, researchers are exploring novel sources and detectors such as QCLs, frequency combs, and metamaterialbased devices, which offer broader spectral coverage and higher resolution [134][135][136][137][138][139] .To expand the bandwidth and achieve sufficient signal strength up to at least 30 THz, further advancements are required, particularly in the development of spintronic 140 and air-based 141 THz generation and detection techniques.Broadband THz spectroscopy can enable researchers to explore a wider range of materials and physical phenomena with higher precision in the future.Another crucial point for THz technology is to increase the signal-to-noise ratio of THz emission.Enhancements in measurement sensitivity are necessary due to the weak absorption of THz radiation by many new materials.Additionally, THz measurements often require trade-offs between bandwidth, signal-to-noise ratio, acquisition speed, and spatial resolution.While increasing the intensity of the THz pulse can improve the signal-to-noise ratio, it can also introduce non-linear and non-equilibrium effects, making the probing process less minimally invasive.Thus, there is a need for novel techniques to improve the signal-to-noise ratio without compromising other measurement parameters.
On the other hand, high efficiency is an important factor in THz spectroscopy, as it directly influences the signal-to-noise ratio and overall performance of the system.Enhancing the efficiency of THz sources and detectors is crucial to improve the sensitivity and reliability of measurements.To this end, researchers are exploring various approaches, including optimizing device structures, enhancing conversion efficiencies, and integrating advanced materials and technologies, to achieve higher efficiency in THz spectroscopy systems [142][143][144][145][146][147][148] .It is worth noting that a recent spintronic heterostructure has demonstrated both a wide THz bandwidth (covering 0.3-30 THz) and high peak electric fields surpassing 1.5 MV/cm.This development holds significant potential in replacing the conventional approach of using LiNbO 3 with tilted pulse-front configuration for generating strong THz fields, offering a more efficient and practical solution 131 .Additionally, a triple-layered chiral structure metasurface was proposed for high-efficiency broadband cross-polarization conversion in THz region, which can reach complete transmission phase coverage (0-2π) at a broadband frequency (0.4-1.0 THz), with 85.7% fractional bandwidth 142 , revealing the potential of integrating metasurfaces and heterogeneous structures for enhanced THz emission efficiency.Similarly, the various metasurfaces and interfaces in heterostructures based on 2D materials probably have great prospects for improving efficiency in the future.Besides, integrability is another key aspect for the practical implementation of THz technology, including trade-offs between various substrate technologies, limitations of device-level metrics, sensitivity to parasitics and packaging issues 130 .Seamless integration of THz sources, detectors, and other components into compact and robust devices would enable portable and field-deployable THz systems for communicating and imaging Integration approaches such as photonic integration, microelectromechanical systems (MEMS), and on-chip technologies are being investigated to achieve compact and integrated THz spectroscopy systems [149][150][151][152] .
In conclusion, the future of THz spectroscopy technology is centered on overcoming key challenges, including achieving broad coverage, high-power operation, high efficiency, and integrability.To address these challenges, there is a need for advancements in THz source and detector technologies, as well as the exploration of innovative materials, heterostructures, and integration approaches, which is driving the development of more versatile THz spectroscopy systems.Furthermore, there is great potential for future 2D materials in THz emission, which encompasses leveraging 2D materials for devices like QCLs to tailor sub-levels for the THz range, along with exploring innovative structures such as 2D materials-photonic crystal heterostructures and Moiré superlattices.These advancements in THz spectroscopy not only expand our understanding of materials and physical phenomena but also open up a wide range of applications in fields such as materials science, biomedical imaging, physics, and communication.

Fig. 1
Fig. 1 THz spectrum interacted with emerging 2D materials.a THz spectrum and its corresponding wavelength and energy.b The Interaction of THz on electrons, spins, and phonons within materials.c vdW heterostructure: based on typical 2D materials.d Schematic illustration of THz emission from 2D materials.e Schematic illustration of THz emission from layered MoS 2 monolayer.f Active THz surface emission in Graphene/SiO 2 /Si heterostructure.g THz emission in WS 2 /MoS 2 heterostructure by optically triggered interfacial current J z .h THz emission induced by Rashba splitting in 2D-phenethylammonium lead iodide (2D-PEPI) crystal.a Reprinted with permission from ref. 1 .b Reprinted with permission from ref.15 .c Reprinted with permission from ref.40 .d Reprinted with permission from ref.52 .e Reprinted with permission from ref.53 .f Reprinted with permission from ref.54 .g Reprinted with permission from ref.56 .Fig h, reprinted with permission from ref.57

Fig. 2
Fig. 2 Schematic of typical experimental setups for THz spectroscopy.a Schematic illustration of typical TES set up. b Schematic illustration of typical THz-TDS set up. c Schematic illustration of typical THz-FDS set up. d Schematic illustration of typical TRTS set up. e THz spectroscopy combines with STM.f Schematic illustration of typical THz-SNOM set up. a Reprinted with permission from ref.58 .b, c Reprinted with permission from ref.65 .d Reprinted with permission from ref.68 .e Reprinted with permission from ref.72 .f Reprinted with permission from ref.78

Fig. 3
Fig. 3 THz emission from non-centrosymmetric 2D monolayers.a Schematic illustration of experimental set up for THz emission from layered 2D monolayer.b Schematic illustration of the resonant optical rectification mechanism for MoSe 2 THz emission.c Time-domain representation of typical THz emissions generated from WSe 2 .d THz emission electric field amplitudes from monolayer WS 2 excited by p-polarization pump beam, with black and red lines representing E x and E y components, respectively.e Time-domain representation of typical THz emission (inset: frequency domain) generated from layered bulk (red solid line) and monolayer (green solid line) MoSe 2 .f Time-domain representation of generated THz emission pulses from layered MoS 2 (red solid curve), graphite (blue dash curve), and InAs crystal (black dash curve).The THz emissions from MoS 2 and graphite are amplified 15 times for better visibility.g Schematic illustration of the THz emission mechanism in FAPbI 3 , showcasing the polaron emission from the photo-Dember effect.h THz time-domain emission spectra of 0.95FA under different pump photon energies (carrier concentration n = 1.25 × 10 18 cm −3 ).a, f Reprinted with permission from ref. 80 .b, e Adapted with permission from ref. 84 .Copyright 2024 American Chemical Society.Fig c, reprinted with permission from ref.82 .d Reprinted with permission from ref.83 .g, h Reprinted with permission from ref.85 THz signal (a.u.)

Fig. 4
Fig. 4 THz emission from photon-drag effect.a Schematic diagram of photon-drag induced THz emission from graphene, pumped by a femtosecond optical pulse.b Comparison between vertical transitions and nonvertical transitions (e.g., from photon-drag).Dashed black contours represent the Fermi surface and green line denotes the energy-momentum-conserving contours.c Spectra in time domain and frequency domain of the emitted THz electric field from monolayer graphene, with incidence angle of opposite signs.d Peak-to-peak amplitude of THz electric field emitted from graphene, as a function of the incidence angle.e THz emission spectra in time-domain obtained from a monolayer graphene sample with a p-polarized pump beam.Red lines, excited from graphene side; blue lines, excited from the substrate side.f Comparison of THz emission from VGG and monolayer graphene with a pumping power of 50 mW.g Dependence of the emitted THz peak-valley value on the QWP angle observed from VGG. h THz-wave magnitudes in frequency domains.Inset: average value of the shaded regions as a function of representative azimuth angles.i Electrically enhanced THz pulses from SnSe 2 crystal.j THz emission signal from the PtSe 2 thin film as a function of time delay under the linear laser pump, compared to the signal from the 0.5-mm-thick ZnTe crystal.k Schematic illustrations of the single-pulse ternary encoding of the polarized THz radiation.a, d Reprinted with permission from ref.86 .b Adapted with permission from ref.87 .c Reprinted with permission from ref.88 .e Reprinted with permission from ref.89 .f Reprinted with permission from ref.90 .g Reprinted with permission from ref.91 .h Reprinted with permission from ref.92 .i Reprinted with permission from ref.93 .j Reprinted with permission from ref.94 .k Reprinted with permission from ref.95

2 Fig. 5
Fig. 5 THz emission from heterostructures.a Gr/SiO 2 /Si and Gr/Si heterostructures' p-polarized THz emission in time domain and frequency domain (inset) with no gate voltage.b Gr/SiO 2 /Si and Gr/Si heterostructures' THz emission at gate voltage -2, 0, 2, and 4 V, respectively.c Graphene/Si Schottky junction p-and s-polarized THz emission amplitude with respect to pump fluence.d WSe 2 /Si heterostructure THz pulse waveform compared with that of Si and monolayer WSe 2 under a femtosecond pulse (2.58 eV and 60 μJ/cm 2 ).e MoS 2 /Si heterostructure THz wave in timedomain from Si metasurface and MoS 2 metasurface in p-polarization state with and without photon pumping (inset is the air reference).f ReSe 2 / MoS 2 and MoS 2 /ReSe 2 heterostructures' THz electric field waveforms (inset: schematic diagram of TES for ReSe 2 -MoS 2 heterostructure) under 800 nm pump excitation (black line for the simulated curve using the photocurrent model).g WS 2 /MoS 2 and MoS 2 /WS 2 heterostructures THz electric field waveforms under 3.1 eV femtosecond pulses (Inset: calibration by a bulk InSb crystal under identical conditions).h THz spectrum analysis of CsPbBr 3 perovskite quantum dots heterostructures under varying laser intensity irradiances.i Co/Bi(n)/Bi 2 Te 3 heterostructures' THz waveform compared with those of ZnTe*0.2,Co, Bi 2 Te 3 , Bi/Bi 2 Te 3 , Co/Bi, Co/Bi 2 Te 3 films.j GaAsSb 0.13 /In 0.141 Ga 0.859 As heterostructure THz waveforms under pump-pulse durations from 140 to 30 fs from red to blue lines.the output THz power P generated by the GL/MoS 2 (k) EDs (inset: GL/b-P) as a function of j. l Temporal evolution of pump-induced THz waveform peak alterations in various samples: graphene/Ag wire grid (black line), CH 3 NH 3 PbI 3 /graphene (red line), CH 3 NH 3 PbI 3 /Ag wire grid/graphene (blue line), and CH 3 NH 3 PbI 3 /Ag wire grid (pink line).The inset in the upper left shows an enlarged view of the peak region.a, b Reprinted with permission from ref.54 .c Reprinted with permission from ref.99 .d Reprinted with permission from ref.98 .e Reprinted with permission from ref.55 .f Adapted with permission from ref.100 , Copyright 2024 American Chemical Society.g Reprinted with permission from ref.56 .h Reprinted with permission from ref.101 .i Reprinted with permission from ref.103 .j Reprinted with permission from ref.104 .k Adapted with permission from ref.97 .l Reprinted with permission from ref.105

Fig. 7
Fig. 7 THz emission from spin-to-charge conversion.a Non-equilibrium electron distribution in the iron layer, excited by a femtosecond laser pulse.b THz emission from Fe/Au and Fe/Ru heterostructures, of which the signal inverts with reversal of the sample magnetization.c THz emission from Co/Pt heterostructure, with the pump is incident from the side of the substrate.d Impact of different NM layer materials on the THz output amplitude of the FM/NM spintronic emitters.e Schematic and the electron population change in the spin injection process from Co to monolayer MoS 2 .f THz emission from Co/MoS 2 heterostructure, compared with pure Co, monolayer MoS 2 and ZnTe crystal.g THz emission from Bi 2 Se 3 /Co heterostructure, pure Co and Bi 2 Se 3 .h THz emission from FGT(4 nm)|TI(8 nm), FGT(4 nm), and TI(8 nm) in frequency domain.i Manipulation of the THz chirality by controlling the magnetic field (Schematic illustration).j Asymmetric spintronic-THz emission from Pt/NiFe and 2D-HMH/NiFe heterostructures (Schematic illustration).k Electric-field-tunable THz spintronic emission under maximum and minimum strained states.l THz emission waves with manipulated amplitudes and polarities.a, b Reprinted with permission from ref.102 .c Reprinted with permission from ref.115 .d Reprinted with permission from ref.116 .e, f Reprinted with permission from ref.117 .g Reprinted with permission from ref.118 .h Reprinted with permission from ref.119 .i Reprinted with permission from ref.120 .j Reprinted with permission from ref.121 .k Reprinted with permission from ref.122 .l Reprinted with permission from ref.123